* ************************************
* ********** Preamble ****************
* ************************************

* Preliminary

clear all
clear matrix
pause off
set more off
prog drop _all

* *********************************
* ****** Tables *******************
* *********************************

use "final_transaction.dta", clear
keep if ln_price_m!=. & treatment_nuclear!=. & property_type_f!=. & rank_deprivation_score!=. & ind_mob_2011!=.

* Table 1: Descriptive statistics in 2010

global variables "price_m property_new property_type_f rank_income_score rank_employment_score distance_nuclear"
global sample "if ln_price_m!=. & year==2010"
su ${variables} ${sample}
global sample "if ln_price_m!=. & year==2010 & treatment_nuclear==1"
su ${variables} ${sample}
global sample "if ln_price_m!=. & year==2010 & treatment_nuclear==0"
su ${variables} ${sample}

* Table 2: Effect of the news shock on property values

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment)
outreg2 using "Table_2", excel label replace keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment lsoa_code)
outreg2 using "Table_2", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_2", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)

* Table 3: The role of neighborhood characteristics--industry composition in 2011

* Columns (1)-(4)

foreach var of varlist *migration_ratio* nomovers_ratio ind_mob_2011 occ_mob_2011 ind_mob_1971 deprived rural unemployed treatmentw_wk children_census share_council_2011 nonwhite_census amenities_public amenities_schools amenities_nature pollution_current {
egen std_`var'=std(`var')
}

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment##c.ind_mob_2011 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_3", excel label replace keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.ind_mob_2011) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment##c.ind_mob_2011 treatment_nuclear##time_treatment##c.occ_mob_2011 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_3", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.ind_mob_2011 1.treatment_nuclear#1.time_treatment#c.occ_mob_2011) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment##c.ind_mob_2011 treatment_nuclear##time_treatment##c.ind_mob_1971 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_3", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.ind_mob_2011 1.treatment_nuclear#1.time_treatment#c.ind_mob_1971) nocons noaster stats(coef se tstat) ctitle(3)
reghdfe ${vary} treatment_nuclear##time_treatment##c.ind_mob_2011 treatment_nuclear##time_treatment##c.deprived treatment_nuclear##time_treatment##c.rural treatment_nuclear##time_treatment##c.unemployed ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_3", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.ind_mob_2011 1.treatment_nuclear#1.time_treatment#c.deprived 1.treatment_nuclear#1.time_treatment#c.rural 1.treatment_nuclear#1.time_treatment#c.unemployed) nocons noaster stats(coef se tstat) ctitle(4)

reghdfe ${vary} treatment_nuclear##time_treatment##c.std_ind_mob_2011 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
reghdfe ${vary} treatment_nuclear##time_treatment##c.std_ind_mob_2011 treatment_nuclear##time_treatment##c.std_occ_mob_2011 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
reghdfe ${vary} treatment_nuclear##time_treatment##c.std_ind_mob_2011 treatment_nuclear##time_treatment##c.std_ind_mob_1971 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
reghdfe ${vary} treatment_nuclear##time_treatment##c.std_ind_mob_2011 treatment_nuclear##time_treatment##c.std_deprived treatment_nuclear##time_treatment##c.std_rural treatment_nuclear##time_treatment##c.std_unemployed ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)

* Column (5)

preserve
use "final.dta", clear
keep lsoa treatment_nuclear treatment_nuclear_pp treatment_powerplant group_nuclear ind_mob_2011 treatmentw_wk rank_deprivation_score residents_census outmigration_ratio_census inmigration_ratio_census migration_ratio_census nomovers_ratio_census
egen tag=tag(lsoa)
keep if tag==1
merge 1:n lsoa using "deprivation_2004_2019.dta"
keep if _merge==3
drop _merge
generate time_treatment=(year>2011)
foreach var of varlist ind_mob_2011 {
egen std_`var'=std(`var')
}
global vary panelrank_deprivation_score
areg ${vary} treatment_nuclear##time_treatment##c.ind_mob_2011 i.group_nuclear##i.time_treatment, absorb(lsoa_code) robust
areg ${vary} treatment_nuclear##time_treatment##c.std_ind_mob_2011 i.group_nuclear##i.time_treatment, absorb(lsoa_code) robust
restore

* Table 4: The role of neighborhood characteristics--moving costs

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment##c.treatmentw_wk ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_4", excel label replace keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.treatmentw_wk) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment##c.children_census ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_4", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.children_census) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment##c.share_council_2011 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_4", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.share_council_2011) nocons noaster stats(coef se tstat) ctitle(3)
reghdfe ${vary} treatment_nuclear##time_treatment##c.migration_ratio ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_4", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.migration_ratio) nocons noaster stats(coef se tstat) ctitle(4)

global vary ln_price_m
reghdfe ${vary} treatment_nuclear##time_treatment##c.std_treatmentw_wk ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
reghdfe ${vary} treatment_nuclear##time_treatment##c.std_children_census ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
reghdfe ${vary} treatment_nuclear##time_treatment##c.std_share_council_2011 ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
reghdfe ${vary} treatment_nuclear##time_treatment##c.std_migration_ratio ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)

* Table 5: Nuclear plants and long-term evolution of neighborhoods (1971–2011)

preserve
use "final_history.dta", clear
reshape long share_labor_oa_ pop_, i(lsoa_code) j(year)
rename share_labor_oa_ share_labor
rename pop_ pop
foreach var of varlist share_labor pop {
egen sd_`var'=sd(`var'), by(year)
egen mean_`var'=mean(`var'), by(year)
generate norm_`var'=(`var'-mean_`var')/sd_`var'
drop sd_`var' mean_`var'
}

global vary norm_share_labor
areg ${vary} treatment_nuclear##i.year i.group_nuclear##i.year, absorb(lsoa_code) robust
outreg2 using "Table_5", excel label replace keep(treatment_nuclear##i.year) nocons noaster stats(coef se tstat) ctitle(1)
global vary norm_pop
areg ${vary} treatment_nuclear##i.year i.group_nuclear##i.year, absorb(lsoa_code) robust
outreg2 using "Table_5", excel label keep(treatment_nuclear##i.year) nocons noaster stats(coef se tstat) ctitle(2)
restore

* *********************************
* ****** Appendix tables **********
* *********************************

* Table A1: Transaction volume

preserve
use "final.dta", clear
global sample "if ln_price!=."
global controls "c.property_type_f##time_treatment c.property_type_d##time_treatment c.property_type_s##time_treatment c.property_new##time_treatment c.property_tenure##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_count
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(robust) absorb(i.month##i.year i.group_nuclear##i.time_treatment)
outreg2 using "Table_A1", replace excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(robust) absorb(i.month##i.year i.group_nuclear##i.time_treatment lsoa_code)
outreg2 using "Table_A1", excel label keep(1.treatment_nuclear#1.time_treatment ) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(robust) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A1", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)
restore

* Table A2: Post-code level treatment

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} postcode_treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A2", replace excel label keep(1.postcode_treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} postcode_treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} i.id_lsoa##i.time_treatment)
outreg2 using "Table_A2", excel label keep(1.postcode_treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} postcode_treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} post_code)
outreg2 using "Table_A2", excel label keep(1.postcode_treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)

* Table A3: Effect of the news shock on property values--quantile regressions

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_f
qreg ${vary} treatment_nuclear##time_treatment i.group_nuclear##i.time_treatment ${controls} ${sample}, quantile(.10)  vce(robust)
outreg2 using "Table_A3", excel label replace keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(.10)
qreg ${vary} treatment_nuclear##time_treatment i.group_nuclear##i.time_treatment ${controls} ${sample}, quantile(.25) vce(robust)
outreg2 using "Table_A3", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(.25)
qreg ${vary} treatment_nuclear##time_treatment i.group_nuclear##i.time_treatment ${controls} ${sample}, quantile(.50) vce(robust)
outreg2 using "Table_A3", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(.50)
qreg ${vary} treatment_nuclear##time_treatment i.group_nuclear##i.time_treatment ${controls} ${sample}, quantile(.75) vce(robust)
outreg2 using "Table_A3", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(.75)
qreg ${vary} treatment_nuclear##time_treatment i.group_nuclear##i.time_treatment ${controls} ${sample}, quantile(.90) vce(robust)
outreg2 using "Table_A3", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(.90)

* Table A4: Pre-treatment differential trends and persistence of the hedonic price response

* Panel A: Placebo check

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment_placebo i.old_new_i##time_treatment_placebo i.tenure_i##time_treatment_placebo c.rank_deprivation_score##time_treatment_placebo"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment_placebo ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment_placebo)
outreg2 using "Table_A4a", excel label replace keep(1.treatment_nuclear#1.time_treatment_placebo) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment_placebo ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment_placebo lsoa_code ${controls})
outreg2 using "Table_A4a", excel label keep(1.treatment_nuclear#1.time_treatment_placebo) nocons noaster stats(coef se tstat) ctitle(2)

* Panel B: Persistence of the hedonic price response

generate time_treatment_duration=0 if year*100+month<201104
replace time_treatment_duration=1 if year*100+month<=201212 & time_treatment_duration==.
replace time_treatment_duration=2 if year*100+month<=201412 & time_treatment_duration==.
global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment_duration i.old_new_i##time_treatment_duration i.tenure_i##time_treatment_duration c.rank_deprivation_score##time_treatment_duration"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment_duration ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment_duration)
outreg2 using "Table_A4b", excel label replace keep(1.treatment_nuclear#1.time_treatment_duration 1.treatment_nuclear#2.time_treatment_duration) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment_duration ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment_duration ${controls} lsoa_code)
outreg2 using "Table_A4b", excel label keep(1.treatment_nuclear#1.time_treatment_duration 1.treatment_nuclear#2.time_treatment_duration) nocons noaster stats(coef se tstat) ctitle(2)
drop time_treatment_duration

* Table A5: Effect of the news shock on property values--large non-nuclear power plants as an alternative control group

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear_pp##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment)
outreg2 using "Table_A5", excel label replace keep(1.treatment_nuclear_pp#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear_pp##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment lsoa_code)
outreg2 using "Table_A5", excel label keep(1.treatment_nuclear_pp#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear_pp##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A5", excel label keep(1.treatment_nuclear_pp#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)

* Table A6: Effect of the news shock on property values--robustness checks with alternative treatment definitions

* Panel A: treated zone

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_10##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6a", excel label replace keep(1.treatment_10#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_15##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6a", excel label keep(1.treatment_15#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_25##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6a", excel label keep(1.treatment_25#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)

* Panel B: 'doughnut' definitions for the treatment

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment ${sample} & (distance_nuclear<20 | distance_nuclear>30), vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6b", excel label replace keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample} & (distance_nuclear<20 | distance_nuclear>40), vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6b", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample} & (distance_nuclear<20 | distance_nuclear>50), vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6b", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)

* Panel C: closer control areas

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment ${sample} & distance_nuclear <= 75 , vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6c", excel label replace keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample} & distance_nuclear <= 50 , vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6c", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample} & distance_nuclear <= 40 , vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A6c", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)

* Table A7: Effect of the news shock on property values--additional controls

foreach var of varlist Green* {
replace `var'=0 if `var'==.
}

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global controls1 "c.Green_vote_2010##time_treatment"
global controls2 "c.rank_income_score##time_treatment c.rank_employment_score##time_treatment c.rank_health_score##time_treatment c.rank_education_score##time_treatment c.rank_crime_score##time_treatment c.rank_barriers_score##time_treatment"
global controls3 "c.house_t_census##time_treatment c.flat_census##time_treatment c.students_census##time_treatment c.retired_census##time_treatment c.managers_census##time_treatment c.professional_census##time_treatment c.high_education_census##time_treatment"
global controls4 "c.distance_shore##time_treatment c.length_rivers##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} ${controls1} lsoa_code)
outreg2 using "Table_A7", excel label replace keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} ${controls2} lsoa_code)
outreg2 using "Table_A7", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} ${controls3} lsoa_code)
outreg2 using "Table_A7", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} ${controls4} lsoa_code)
outreg2 using "Table_A7", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(4)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year##i.group_nuclear ${controls} lsoa_code)
outreg2 using "Table_A7", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(5)

* Table A8: Effect of the news shock on property values--Nationwide data

preserve
use "final_nationwide.dta", clear
cap egen id_lsoa = group(lsoa_code)
global sample "if lnprice!=. & rank_deprivation_score!=. & year>=2007"
keep ${sample}
global controls "c.bedrooms##time_treatment c.floorsz##time_treatment c.bathroom##time_treatment c.rank_deprivation_score##time_treatment"
global cluster_list "id_lsoa month"
global vary lnprice
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, absorb(i.month##i.year i.group_nuclear##i.time_treatment lsoa_code) vce(cluster $cluster_list)
outreg2 using "Table_A8", excel label replace keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment ${sample}, absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code) vce(cluster $cluster_list)
outreg2 using "Table_A8", excel label keep(1.treatment_nuclear#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_waste##time_treatment ${sample}, absorb(i.month##i.year i.group_waste##i.time_treatment lsoa_code) vce(cluster $cluster_list)
outreg2 using "Table_A8", excel label keep(1.treatment_waste#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)
reghdfe ${vary} treatment_waste##time_treatment ${sample}, absorb(i.month##i.year i.group_waste##i.time_treatment ${controls} lsoa_code) vce(cluster $cluster_list)
outreg2 using "Table_A8", excel label keep(1.treatment_waste#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(4)
restore

* Table A9: Effect of the news shock on property values--nuclear waste sites

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_waste##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_waste##i.time_treatment)
outreg2 using "Table_A9", excel label replace keep(1.treatment_waste#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_waste##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_waste##i.time_treatment lsoa_code)
outreg2 using "Table_A9", excel label keep(1.treatment_waste#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_waste##time_treatment ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_waste##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A9", excel label keep(1.treatment_waste#1.time_treatment) nocons noaster stats(coef se tstat) ctitle(3)

* Table A10: Effect of the news shock on property values--the role of plant characteristics

egen active_nb_census_20=sum(active_nb_census*(treatment_nuclear)), by(group_nuclear year month)
generate nb_employees=750/active_nb_census_20 

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment##c.connection ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A10", excel label replace keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear##1.time_treatment#c.connection) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment##c.expclosure ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A10", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear##1.time_treatment#c.expclosure) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment##c.numberaccidents ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A10", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear##1.time_treatment#c.numberaccidents) nocons noaster stats(coef se tstat) ctitle(3)
reghdfe ${vary} treatment_nuclear##time_treatment##c.packvol_waste ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_waste##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A10", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.packvol_waste) nocons noaster stats(coef se tstat) ctitle(4)
reghdfe ${vary} treatment_nuclear##time_treatment##c.nb_employees ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_waste##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A10", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.nb_employees) nocons noaster stats(coef se tstat) ctitle(5)

* Table A11: The role of neighborhood characteristics--amenities and geography

* Panel A: amenities

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment##c.amenities_public treatment_nuclear##time_treatment##c.amenities_schools ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A11a", excel label replace keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.amenities_public 1.treatment_nuclear#1.time_treatment#c.amenities_schools) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment##c.amenities_nature ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A11a", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.amenities_nature) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment##c.pollution_current ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A11a", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.pollution_current) nocons noaster stats(coef se tstat) ctitle(3)

* Panel B: geography

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
reghdfe ${vary} treatment_nuclear##time_treatment##c.mean_elevation ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A11b", excel label replace keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.mean_elevation) nocons noaster stats(coef se tstat) ctitle(1)
reghdfe ${vary} treatment_nuclear##time_treatment##c.orientation_nuclear ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A11b", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.orientation_nuclear) nocons noaster stats(coef se tstat) ctitle(2)
reghdfe ${vary} treatment_nuclear##time_treatment##c.length_rivers treatment_nuclear##time_treatment##c.distance_shore ${sample}, vce(cluster $cluster_list) absorb(i.month##i.year i.group_nuclear##i.time_treatment ${controls} lsoa_code)
outreg2 using "Table_A11b", excel label keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.length_rivers 1.treatment_nuclear#1.time_treatment#c.distance_shore) nocons noaster stats(coef se tstat) ctitle(3)

* Table A12: The role of neighborhood characteristics--interaction terms

global sample "if ln_price_m!=."
global controls "i.property_type_i##time_treatment i.old_new_i##time_treatment i.tenure_i##time_treatment c.rank_deprivation_score##time_treatment"
global vary ln_price_m
global cluster_list "id_lsoa month"
global controls6 "treatment_nuclear##time_treatment##c.pollution_current treatment_nuclear##time_treatment##c.amenities_nature treatment_nuclear##time_treatment##c.amenities_public treatment_nuclear##time_treatment##c.treatmentw_wk treatment_nuclear##time_treatment##c.share_council_2011 treatment_nuclear##time_treatment##c.migration_ratio_census treatment_nuclear##time_treatment##c.children_census treatment_nuclear##time_treatment##c.deprived treatment_nuclear##time_treatment##c.unemployed treatment_nuclear##time_treatment##c.rural"
reghdfe ${vary} treatment_nuclear##time_treatment##c.ind_mob_2011 i.month##i.year i.group_nuclear##i.time_treatment ${controls} ${controls6} ${sample}, vce(cluster $cluster_list) absorb(lsoa_code)
outreg2 using "Table_A12", excel label replace keep(1.treatment_nuclear#1.time_treatment 1.treatment_nuclear#1.time_treatment#c.ind_mob_2011 1.treatment_nuclear#1.time_treatment#c.pollution_current 1.treatment_nuclear#1.time_treatment#c.amenities_nature 1.treatment_nuclear#1.time_treatment#c.amenities_public 1.treatment_nuclear#1.time_treatment#c.treatmentw_wk 1.treatment_nuclear#1.time_treatment#c.share_council_2011 1.treatment_nuclear#1.time_treatment#c.children_census 1.treatment_nuclear#1.time_treatment#c.migration_ratio_census 1.treatment_nuclear#1.time_treatment#c.deprived 1.treatment_nuclear#1.time_treatment#c.unemployed 1.treatment_nuclear#1.time_treatment#c.rural) nocons noaster stats(coef se tstat) ctitle(1)

